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Abstract 

In the paper the interfering resonances parameters determination ambiguity is considered. 
It is shown that there are two solutions for two fixed width resonances. Analytical relation 
between different solutions is derived. Numeric experiments for fixed width three and four 
resonances, and for model energy-dependant width two resonances confirm ambiguity of the 
resonances parameters determination. 

I. INTRODUCTION 

One of the typical tasks during experimental data processing is determination of the 
parameters of several resonances from the experimental cross-section measurements taking 
into account their interference with arbitrary phases. Often it occurs that for the resonances 
with arbitrary phases several almost equally good solutions can be found. In the paper an 
attempt to analyze this problem is performed using simple examples. 

The preliminary version of this paper has been published as a preprint of Budker Institute 
of Nuclear Physics [l| . 

II. TWO RESONANCES INTERFERENCE WITH SIMPLE 
NON-RELATIVISTIC BREIT-WIGNER AMPLITUDE 

Let us consider a model cross 
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where A, B are some complex numbers (coupling constants), mi, m2, Fi, r2 are real numbers 
(masses and widths of resonances) . 

First let us derive the conditions which lead to the identical cross section as a function 
of energy with different set of parameters. 

Two identical continuous functions should have identical Fourrier images. 

Fo. the function of .nt.est the Fou.*. i.age Q is easily calculated: 
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In order that function 
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be equal to c^E) at every point E, evidently the following equalities should be valid 
mix = mi, F13. = Fi, 
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Apparently the resonance masses should be ordered here otherwise additional trivial solu- 
tions would appear due to parameter sets exchange. 

For the amplitudes Ax, Bx we have four equations with four variables (separate equations 
for the real and imaginary parts). Because the equations are non-linear there could be more 
than one solution. 

Since only amplitude absolute value squared has a physical sence there is a freedom in 
absolute phases with definite relative phase value. So we can take, for example, that Ax is 
a real number and Bx defines their relative phase, or equivalent: 

Ax I Ax I 6 , Bx I Bx I 6 , 
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If the latter definition is admitted then 



(5) 



where a, b, ip, a^, b^, ipx are real numbers. 

Now one can write down system of equations: 



^ O'xOx-] Q 



(ri+r2)+{mi-m2)sin{2i/)^)] 



a 



a6-[i2^2^(ri+r2)+(mi-m2)sm(2i/)) 



"°°'''^'"' (ri+r2)-(m2-mi)sin(2i/^^)] _ 



f)2 axbx-^ 2 

^ ^;mi-rn2)V(^)' 

_ ab.[i^2Sm(ri+r2)-(m2-mi)sin(2^)] 

ax bx ■ [(mi -m2 ) cos(2fe )- ' (Fi +r2 )] _ 

_ afe-[(mi-m2)cos(2V')-5iii^(ri+r2)] 

~ (,ni-m2)V(^y 

axfex- [(m2 -mi ) cos(2fe )+ ' (Fi H-Fa )] _ 

(,ni-m2)V(^y ~ 

_ a6-[(m2-mi)cos(2V')+5iil^(ri+r2)] 



(6) 



(mi-m2)^+(^)' 

Evidently the last two equations are identical, so for three unknown variables ax,bx,ipx we 
have three independent equations. 

Trivial solution: = a, b^ = b, ipx = ^- Let us check whether there are some other 
solutions. First exclude ip^- One equation without ip^ can be derived by subtracting the 
second equation in (JH]) from the first one: 



at b 
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Let us introduce a new variable 



al - a? bl- b'^ 
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Now for two variables ip^ and y we have two equations: 



cos(2Vj^) 



(ri + T2) + (mi - 7722) sin(2?/;j 



(Fi + Ts) + (mi - ma) sin(2^A) 
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(mi - ma) cos(2^.) - (Fi + T2) 



ah ■ 



(mi - ma) cos(2^A) - ^ (Fi + Ta) 



Linear equation for tg{2il)x) can be obtained dividing the first equation in ([9]) on the 
second one: 



^^4^+(mi-m2)tg(2V>^) 
(mi-m2)-M|M(ri+r2) 



a6- 



^(ri+r2)+(mi-m2)sin{2^)]-y[(mi-m2)2 + (£4ii)'] 
a6.[(mi-m2)cos(2V')-2iii^(ri+r2)] 
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The solution is amazingly simple: 



a b sin(2'?/;) + (ma — mi) y 



(11) 



a6cos(2^)-|(ri + ra)' 
Now we can check whether it is an actual solution of the system (Q. Let us check the 
following values: 



Si 



ax ■ (Fi +r2 )+ (mi -m2 ) sm(2V>x )] 



ab.[£2£m(ri+r2)+{mi-m2)sin{2V0]-y[(mi-m2)'+(^4^) ] 



_ a,6,-[(mi-m.2)cos(2V:r)-^i^i^^(ri+r2)] 
^ ~ afe-[(mi-m2)cos(2V')-2i2^(ri+r2)] ' 



(12) 



For an actual solution there should be S'l = 1, ^a = 1. After substitution tp^ we obtain: 

a^b^ cos(2V'a;) 



iSi — S2 



a6cos(2^)-|-(ri + ra) 



From S'l a = 1: 



cos(2?/';j 
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a6 cos(2V') - f ■ (Fi + Fa) a,6, ^(a^ + yFi) (62 + yFa) 
and from expression flTT]) we get: 

cos(2-0x) 

a6cos(2i/i)-|-(ri+r2) ~ 

_ 1 

ya2b2+»/^-[(mi-m2)^ + ( '"l3^^ )^]+2aby-[(m2-mi)sin(2t/>)- '^i+'^2 cos(2^/))] 
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So the equation for y: 



(16) 



+ 2a 6 ■ [(m2 - mi) sin(2^/;) - cos(2V^)] - a^Fs - h'^Vi ] = 
This equation has two solutions for y: 



y = 



a, = a, K = b, tg(2^,.) = tg(2^), 8^ = 82 = 

It is obvious that the values sin(2?/;j.) and cos{2tlJx) must match exactly with sin(2?/^) and 
cos{2ip), correspondingly. 



y 



a2r2+62ri+2ab-[22£m(ri+r2)+(mi-m2)sm(2V') 



\/a?~+yI\, = \/lP~+yI^, 
a6cos(2V')-|(ri + r2) 



(^x bx 

absm{2il)) + (m2 — mi)y 
^x bx 



(17) 
(18) 
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If one formally substitutes this solution, then two cross section curves become identical. 
However in order to this solution be admittable some conditions must be satisfied: 



I. y>-f^^al>o 



2. y>-f^^bl>0, 

3. |cos(2V'a:)| < 1, 

4. \sm{24j,)\ < 1. 



Let us check 



0^ = 0^ + Tiy 



fe2r2+a2.[(mi-m2)2 + (i^)']+2aferi.[^2;5m(ri+r2)+(mi-m2)sm(2V')] 



2). 



(20) 



(mi-m2) + 
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ft^r? cos2 ( 2^+arctg ^,'"i+^2 

i V ^ 2{mi ~m 



2), 



(mi-m2) + 



2:(El-l2' 



> 0. 
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Similarly the condition b'^ > is checked. The last two conditions are easily confirmed by 
the check that derived by their own formulae sm[2iljx) and cos(2-?/^a.) satisfy the Pythagorean 
theorem cos'^{2ipx) + sm'^{2tpx) = 1- 

So we have found non-trivial solution that means: any pair of resonances can be replaced 
by another pair of resonances with the same masses and widths but different amplitudes 
and phases so that the resonance cross section does not change. 

Using these formulae one can get the solution for the case of interference of Breit-Wigner 
amplitude with complex constant, substituting 



and setting r2 ^ oo. Corresponding cross section reads 



ae 



+ ce 



-iip 



(21) 



Finally we get 



- lim "''^2 + (c-iii)^ri+a-c-j-r2-[22£m(ri+r2)+(mi-m2)sin{2i/;)] 

= 2iaccos2ip — c^Fi, 



(22) 



a 



^a? + 2iacTi cos 2^ - c^F^, 
ig2i)x = "t'^-v 

Here the variable y, which was real in previous consideration became complex as well as 
variable ax- If reassemble imaginary and real parts of these variables or solve this problem 
from the beginning with cross section ( !2T!) . then one gets 



(23) 



'a^ + 2acFi sin 2il) + c^Ff = ^ (cFi + a sin 211))' + cos2 2^^, 

sin 27/;, = cos 27/;, = ^i^. 

Let us look at the numerical example of uo and mesons interference in the channel 
e^e~ n^n^Ti^. Approximate values of resonance parameters: mi = m^^ = 782.6, F^^ = 8.4, 
a = i^m2 = m^ = 1019.4, F^ = 4.46, b = 0.1, ip = = -77.5°. 

Substituting to the formulae one gets 



y = 0.00042; a, = 1.0018; 6, = 0.109; t/jx = 74.4°. 



(24) 
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So we get quite different phase while amphtudes changed just a httle bit. 

At this point the most urgent question is: whether this ambiguity is a unique property of 
just this simple resonance description, or in a more sophisticated and realistic parameteriza- 
tion of resonance cross section such resonance phase ambiguity will take place? The matter 
is that for actual experimental data processing the much more complicated resonance cross 
section formulae are used. 



III. RELATIVISTIC BREIT-WIGNER RESONANCE AMPLITUDE 



A little more complicated variant of Breit-Wigner formula (relativistic) is: 



aiE) 



A B ^ 



+ 



(25) 



s — m{ + iTimi ' s — m| + iT2Tn2 
where s = E"^, evidently has the same property, just some redefinition of variables is neces- 
sary. 

It is also evident that additional general factor, even strongly dependent on energy, does 
not change the solution. 

However in the most accurate variant of this formula instead of constants Fi, r2 there 
are used more complicated expressions containing phase space of final states, transition to 
which are probable for these resonances. 



IV. ENERGY DEPENDENCE OF RESONANCE WIDTH 



For multiparticle final states, such as vr+vr 7r°, there are no simple formulae for the final 
state phase space, so for our exercise let us choose a simple model formula: 



s - Afi"^ 



(26) 



where the effective mass fi defines the reaction threshold, so at s < cross section becomes 
equal zero. New cross section can be rewritten as follows 



a{s) 



( s—4fi 



, 3 
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s—mf+iTimi 
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For this function it is hard to make a Fourrier transform in order to apply the same trick 
we have done for the approximate resonance curve. Substitute 



s = 4:fi^ + p^, p 



V)4 



pi = {ml - V)4 ^ p2 = {m'i - V 
Cross section dependence on the new variable p is the following 



(28) 



■X 



X 



P2 



(29) 



The task of Fourrier transform of this function is already not so hard if we can find all 
its irregular points. Just let us simplify the function first so as the general factor does not 
change the problem solution. 

Thus investigated function of p is 



/(P) 



+ 



B 



p-l+4/i2-m|+ir2m2- 



,P2 



(30) 



16 irregular points are determined by the equations 

3 



p + 4p — m • ± iTjirij ■ 



0, J = 1,2 



Rewrite this equation in the form 



i?iP^ + i?2 = 0, 



(31) 



(32) 



where Ri = ±i--r-, --^z -r- ■■-j rj 



4p 



Pt 



If we solve this equation in a standard way, then the solution can be excessive complicated. 
Let us try to apply the idea of Ferrari solution directly to the equation fl32|) . 



p^ + Rip^ + R2= (p^ + Aip + As) ■ (p^ + Asp + A4) 
This equality is valid for any p if 

Ai + A3 = Ri 
A2 + A4 + A1A3 = 
A2A3 + A1A4 = 
A2A4 = R2 



(33) 



(34) 
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Exclude Ai, A3: 

Ai = -^i?i, A3 = -^i?i. (35) 



Now for the variables A2 and A4 we have two equations: 

(A2 + A4) ■ (Ai + Al - 2R2) - RIR2 = 0, 

A2A4 = R2- 

Substituting A2 + A4 = y, obtain 



Ai 



A 



2 



A4 



Substituting Ri, R2, we get 



lit I / 



(36) 



(37) 



and ?/ must satisfy the equation 

- 4R2y - R\R2 = 0. (38) 

Applying Cardano solution 

y = a + p ^ + + {a + p) ■ {Sap - AR2) - RIR2 = 0. (39) 
Setting P = we get 

,6 d2d -i , ^ 



- RiR2a' + ( ^ ) =0, (40) 



(41) 



(42) 



It can be easily confirmed that a > 0, /? < 0, ?/ > 0, A2 > 0, A4 < are the real numbers, 
Ai, A3 are complex numbers. 
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Two roots are determined by the equation 

p2 + Aip + A2 = 0, (43) 

and two more roots by the equation 

p2 + AsP + A4 = 0. (44) 

None of these roots can be a real number that is evident from the initial form of quatric 
equation: p = is not root, and for any real p 7^ polynomial has non-zero imaginary part 
for non-zero width Tj > and mass rrij > 0. In order to calculate integrals with infinite 
limits by residue method we are interested to know the sign of imaginary part of roots. 
Introduce notation for all roots. 



~* o A 2 , ..4 




(45) 





— A2I: 




— A2I 




— A4I, 


''^31 


— A4I 



~ A31 

Z41 2- 

One can see that ^{zu) = ^{z2i) = 0, '^{zu) > 0, $5(221) < 0, '^{z^) < 0, '^{z^i) < 0. 
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Four roots correspond better to the equation with index j = 2: 

X _ ^- y2+\/yi+4p| Tama 



A22 = I (z/2 + VVI + 4^2) 
A _ „- j;2-A/j;2+4P2 Tama 

A42 = i (z/2 - Vl/I + 4p^ 




+ (ir 



^12 


^ ^22; 


''^12 


— A22: 


-^^32 


— A42; 




— A42 



~ A12 

2^12 2~ 

~ A12 

^22 2~ 

~ A32 

Z32 2~ 

~ A32 

2:42 2~ 

Here ^{z^) = 3^(^22) = 0, Q{zi2) > 0, Q{z22) < 0, 53(^32) < 0, ^(^42) < 0. 
Now the function /(p) can be presented in the form 

ftp) = ( A \ B \ ^ 

•' \{P-Zll){p-Z2l){p-Z3i)(p-Z4l) (p-Zi2)(p-222)(p-232)(p-242) / 

X I , — ^ — 4^ — ^ — , + 



(p-^n) (p-4i) (p-^3*i) (P-^1*2) (p-42) (p-42) (p-^12) 



(46) 



(47) 



where symbol * designates complex conjugate number. Fourrier image F{t) = J /(p)e**'' dp 

—00 

of the real function has a property that F{—t) = F*{t), so it is enough to calculate Fourrier 
transform only for positive value of t, that is determined by the sum of residues on the 
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irregular points above the abscissa axis, that is 2:11, Z21, z^^, z*^-^, Z12, Z22, -^32, zl. 



F{t) = 2mx 



{Zll-Z2l){zil-Z3l){zil-Z4l) \^ 2211 (211-221) (211 (211 

+ m 

(^11-42) (211-42) (^11-42) (211-^42) 



+ 



(4i-2n)(4i-4i)(4i-4i) y(^2i-2ii)24i(4i-23i)(4i-24i) 

B 



(4i -212) (41-222) (41-232) (41-242) 



+ 



(2^1 -2n ) (2^1 -2J1 ) (2^1 -2j^ ) I (4l -Zll ) (4l -Z21 ) (4i -Z-il ) (2*^ -241 ) 



B 



(4l-Zl2)(z5i-Z22)(4i-^32) (41-242) 



_| A*e'*^4i _ / A |_ 

( 2| 1 - Zn ) ( 2|i - 2|i ) ( 2|i - 2 J J I (2I1 -211 ) (2I1 -221 ) (zli -Z31 ) (4l -241 ) 



B 



(4l "212 ) (z|i -222 )(zJi -232) (41-242) 



+ (48) 



_l ge'*^i2 _ / A^ I 

(Z12-222)(212-232)(212-Z42) I (^Zl2-Z*^)[zi2-Z*^)[zi2-Z^i) (zi2-2| J 



B* 



2212(212-2^2) (212-42) (212-2I2) ^ 

_B*e'*42 / A 



(2^2-2^2) (222-232) (252-2J2) V (222-211) (2^2-221) (2J2 -231) (2J2 -241) 



B 



~'~ (2J2 -212 ) 22*2 (2J2 -232) (2J2 -242 ) _ 



(42-2^2) (2^2-222) (2I2-2J2) \ (232-211) (2I2-221) (2I2 -231) (2^2-241) 
_^ B j _^ 

(2*2-212) (2*2-222) (2^2-232) (2*2-242) y 

_| B*e'*42 _ / A ^ 

(2^2-2^2) (2I2-232) (2^2-212) \ (242-211) (2|2-22l)(2|2-23l)(2|2-24l) 

I B 



(2I2 -212 ) (2I2 -222 ) (2I2 -232 ) (2I2 -242 ) 

Because the coordinates of irregular points determine the functional dependence of Four- 
rier image on the parameter t, so they must be unchanged for any solution. This condition 
produces the equations: 

+ Afil - m% ± tTj^nij, ( = + V - ^? ± ^^,^3 f-) (49) 

V Pjx / \Pj / 
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for j = 1,2. Thus we get the system of equations: 

Tixjuix _ rimi 
T'2xrn2x _ r2m2 



(50) 



where pi = {m\ — Ajj?)^^^ , p2 = {ml — Ajj?)^^^. For five variables fix,^ix,^2x,mix,m2x we 
got the system of the four equations. However the factor in the original function de- 
termined the threshold behaviour, and from p = {s — Ap?)^^^ is evident that p must not 
change, so 

Px = p (51) 

and we have four equations for the four variables. Taking into account that mj^. > 0, we 
derive mjx = mj and can immediately conclude that pjx = pj, and Tj^ = Tj. 
Let us try to simplify the denominators in the formula (HHl) . Apparently 

2^11 (2:11 — 2:21) (2^11 ~ 2:3;^) {zii — z^-^^) = 

/ \ 3 

= + 4/1^ -mj- iVimi I ^ 
= zf^ + 4/1^ -m\- iVimi 

zf^ + 4/^2 -ml + iTimi (^^y = -2iTimi 



(52) 



Similarly 



.z 



21 - zii) 2^21 iz*i - ^31) i4i - ^41) = 2irimi ( — 

Pi 



'31 



^11) (2^31 - 2:21) (^31 - 231) (2^31 - ^41) = 2iVimi ( — 

Pi 



[zli - zii) {zl^ - Z21) (4i - ^31) (4i - ^41) = 2iTimi 
2zi2 {zi2 - Z22) (^12 - 232) (212 - ZI2) = -2iT2m2 



£41 
Pi 



Z12 

P2 



.Z 



22 ~ ^12) 2^22 (^22 ~ ^32) (^22 ~ ^42) — 2ir2m2 ( — 

P2 



^32 ~ ^12) (2^32 ~ ^22) (^32 — 2^32) (2^32 — 242) — 2iT2m2 ( — 

P2 



(53) 
(54) 
(55) 
(56) 
(57) 
(58) 
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(^2 - ^12) (^2 - 2:22) (^42 - ^32) (^2 " ^42) = 2zr2m2 (^^^ (^9) 



Eight other denominators are not so compact: 

(zii — (^11 ^22) (^11 ~ ^32) (^11 ~ ^42) 



-zr2m2(^) -zrimi(^ 



(^21 ~ ^12) (^21 ~ ^22) (^21 ~ 2:32) (2:21 ~ ^42) 

3 / * \ 3 



^^31 ~ ^12j 12:3^ - 222^ 1^31 - 2:32; ^^31 - 2;42j 

3 / * \ 3 



^r2m2 (^^) + iTirrii (^^^ -ml + ml 



[Z^l — Z12) [Z^i — Z22) — Z32) [^4^1 — Z42) — 

= iT2m2 (^^) + iVinii j -ml + ml 

{Zi2 - Zli) {Zi2 - 221) (2^12 - Zli) {Zi2 - 241) = 

= -2r2m2 (^)' - iT^im i^^y -ml + ml 

(^22 ~ ^11) (^22 ~ ^21) (^22 ~ ^31) (^22 ~ ^41) = 

= ir2m2(^^) +iTimi(^^^ -ml + ml 



{ZI2 - Zii) (Z32 - Z21) (^32 - 2^31) (^2 - ^41) 

3 / . \ 3 



iV2m2 + iT\mi [^^^ — ml + ml 



(^42 — 211) (^42 — 2:21) (2:42 — 231) (^42 ~ ^41) — 

= ir2m2 (^^) + iTimi (^^^ -ml + ml 
Here we can shorten the expressions using notations: 



(60) 



(61) 



(62) 



(63) 



(64) 



(65) 



(66) 



(67) 



Timi Tama 2 2 2 /«qx 
G = ^^ + ^^, D^ = m2-m^. (68) 

Pi P2 

For the rest of free parameters Ar^ = axC^^"^, = h^e'^^^ there is a system of eight 
complex equations: 



ae 



-iiP I ae'^ _|_ bej''^ 



2jrimi ( 



P2 



22 "^"^ s^r^maf^'' 



11/. be~'^ 



2iT2m2( 
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iGzll+Dl, 



Let us introduce new variables: 



5e# 



ah 







2iT 2 ,112 









(76) 



+ 



V 



(77) 



Of course the trivial solution of the system of equations is valid: ipx = = a, bx = 

b, ^ = 1, Ux = y. 

The new set of unknown variables ipx^ ^x-, ^- Old variables are connected with new ones 



with the following relations: = a/ ab^ Vx , b^ 



The system of equations with new unknowns ^, Vx is: 



~7 1 7 

^ 2irimi(^ 



-2il/)a; 



-2il/l;c 



-211/" 



PI 



-2ii/> 



+ 



2irimi 



21 
PI 

3_n2 



2jrimi 



^31 
PI 

2ii/' iGz't'^ — Hp- 



-2ii/> 



2irimi ^ 
2iV2m.2[^ 



-,2iipx 



+ 



2iT2m2[^-^ 
iGz*i+Dl 



2ir2m2 



\-yx 



\--x 
\ ]_ 

Cl' Vx 



\ \_ 

iv Vx 

\ \_ 

i,V Vx 

\ \_ 

Vx 



(78) 



Subtracting the third equation from the second one, we get 



2iTimi 



2iTimi 



\ 



I 



(79) 



which can be only valid if one of the multipliers equals zero. So one of the solution is 



(80) 
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It can be checked easily that the first multipher cannot be equal to zero for any resonance 
parameters. If it were so then 

^G--f = ^G--^, (81) 

^21 ^31 

and this is equivalent to the equality Z21 = Z31. Similarly using the equation pairs "second 
- fourth", "third - fourth", and supposing 7^ |, one gets 221 = ^31 = -2^41- So the three 
roots of quatric equation are equal to the same value. 

p4 + 4/i2 -ml + ilimi (^-^y = {p- zu) [p - Z2if = 

= p^-(z,, + 3z2i)p'+ (82) 

+ 82:21 (2:11 + 2:21) - 4l (3^11 + 2:21) p + 2ii2^i. 

So as none of the roots is equal to zero, then comparing to representations of the same 
equation one concludes that 

Zii + Z21 = 0, 

(83) 

3zu + Z21 = 0, 

and hence zu = Z21 = 0, which is impossible. So our guess 7^ | is invalid and the solution 
of ([79D is 

= J. (84) 
Similar analysis of equations 6-8 in fl78|) gives 



= ^u. (85) 

So one can immediately conclude that 

e = i, '^. = '^- (86) 

Now from any equation from the system flTHl) one derives 

^2i^. = e^^^^ ^ = ^. (87) 

As a result we conclude that for the case of energy dependent resonance width the dege- 
naration disappears and there is the only set of resonance parameters presenting the given 
energy dependence of cross section. 
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V. NUMERIC EXPERIMENTS 

For additional check of these conclusions it is useful to carry out numeric experiments 
simulating some actual experiment in high energy physics. Such experiments can show the 
role of experimental statistics. 

In all cases we should suppose some true process cross section function <j{E), where 
E is the energy of colliding beams, then for the finite number of points we generate 
experimental number of events, according to the Poisson probability distribution (integrated 
luminosity at every point is equal to the same value L). 

In order to get the parameters of resonances, consisting the given model of process, we 
minimize the likelihood function as follows 

(88) 



£ = ^ 2 ( pfe - nfc + rifc In — j 
fc=i ^ 



which is equal to the doubled logarithmic likelihood function with opposite sign. Here 
Pk = <^{Ek) ■ L. For greater statistics this function limits to value, that can be used for 
the check of statistical confidence level. Minimization will be performed with well-known 
MINUIT package fl. 

For presentation of experimental cross section on the plots, the experimental number of 
events will be ascribed asymmetric statistical errors: 

An5+^ = V^^~Tl, Ant^ = V^. (89) 
A. Approximate expression for the resonance amplitude 



a{E) = 

with "true" values of parameters 



E — nil + E — 1712 + 



(90) 



mi = 782.6, Fi = 8.4, a = 1, 

(91) 

7712 = 1019.4, F2 = 4.5, 6 = 0.3, = 155°. 

This function is equal to a{mi) ^ ^ = 0.057, 0"(m2) ~ = 0.018. In order to an 
accuracy at the resonance maxima to be at least at the 5% level let us appoint L = 2 ■ 10^. 
Fig.[T] shows the "experimental" set of points. 
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800 900 1000 E 800 900 1000 E 

Figure 1: Result of fit to "experimental" points. Cross section model is ([90]) . x^/^D = 64.3/(74—7). 




\)/2^, degrees 

Figure 2: Plot of the likelihood function on the phase of the second resonance amplitude ■023; 

Plot of likelihood function on the phase of second resonance amplitude ip2x is presented 
in Fig[2l Two equivalent minima with x^/^-D = 56.703/(74 — 7) are obtained at ip2x = 
— 157.14° and ip2x = 157.51°. Although the separating maximum is not high (x^ = 56.907 at 
'iIj2x = ±180°), these are the different solutions. Values of all parameters at these minimum 
points are cited at the Table[Tl 

This numeric experiment confirmed the analytical conclusion — fitting the data with the 
approximate resonance formula produce the ambiguity of resonance phases and amplitudes. 
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Table I: Parameters of resonances at the minimum points of likelihood function. 





mix 


Tlx 


bx 


'm2x 


^2x 


r2x 


x' 


1.0125 


782.62 


8.5371 


0.30284 


1019.4 


4.5554 


157.815 


56.703 


1.0167 


782.62 


8.5371 


0.31033 


1019.4 


4.5554 


-157.065 


56.703 



Table II: Resonances parameters at the two minimum points of likelihood function. 





mix 


Tlx 


bx 


m2x 


^2x 


r2x 




1.0020 


782.60 


8.4116 


0.30007 


1019.4 


4.5093 


155.283 


74.58624 


1.0071 


782.60 


8.4116 


0.30708 


1019.4 


4.5093 


-154.807 


74.58624 



B. Relativistic form of resonance amplitude 

Let us consider more accurate dependence of the resonance amplitude on energy: 



a{E) 



+ 



E"^ — ml + iTirrii E"^ — m\ + iT2m2 
with the same "true" parameters values 



(92) 



mi = 782.6, Ti = 8.4, a = 1, 

m2 = 1019.4, 12 = 4.5, b = 0.3, = 155= 



(93) 



where the factors 2mj are introduced in order to keep the valeus of amplitudes at the 
resonance masses, and general factor imitates the threshold behaviour of cross section with 
^ ^ 3rnjr_ 140+140+135 ^ 2O8. In Order that minima of likelihood function to be more 
demonstrative let us increase the integrated luminosity to the value of L = 10^. 

In Fig. [3] the set of points and optimal cross section of the type (1921) are presented. 

The plot of the likelihood function on the second resonance phase ip2x is presented in 
Fig.l 

As we expected, we get two equivalent minima again. Parameters values at the minimum 
points of C are shown in TablellTl 
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800 900 1000 E 800 900 1000 E 

Figure 3: Result of fit of "experimental" points. Cross section model is described with the for- 
mula (IMD. X^luD = 74.6/(74 - 7). 




\|/2„ degrees 



Figure 4: Plot of likelihood function on the second resonance phase '02a; 
C. Resonance width dependent on energy 

Finally let us consider the case where the degeneracy is expected to disappear and only 
one global minimum will be found: 
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800 900 1000 E 800 900 1000 E 



Figure 5: Result of fit of "experimental" points. Cross section is described by the formula (I94p . 
xV^Z) = 74.8/(74 -7). 

Table III: Resonance parameters at the two points of minimum of likelihood function. 



ax 


mix 


Tlx 


bx 


m2x 


^2x 


r2x 


x' 


1.0020 


782.60 


8.4117 


0.30007 


1019.4 


4.5093 


155.289 


74.82254 


1.0071 


782.60 


8.4117 


0.30699 


1019.4 


4.5093 


-155.529 


74.81663 



X 



2mia 



E'^—mf+iTimi' 



2m2be*^ 



E'^—m%+iV2m2- 



(94) 



Fig. [5] shows the set of energy points and optimal cross section of the type ( l94l) . 
Likelihood function plot on the phase of the second resonance iIj2x is shown in Fig. [61 
Indeed, the minimum values are not equal now (see Table HTTll . however the difference 

is very small and, what was totally unexpected, the "better minimum" corresponds to the 

"wrong" minimum. 

Let us look what will change if the experimental statistics will increase by factor of 100 
(L = 10®). In Table [IV]the parameters of resonances are shown for two minimum points 
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\|/2„ degrees 

Figure 6: Likelihood function plot on the phase of the second resonance ip2x 
Table IV: Resonance parameters at the two minima of likelihood function. 





mix 


Tlx 


bx 


m2x 


^2x 


r2x 




0.99982 


782.60 


8.3977 


0.30015 


1019.4 


4.5025 


154.985 


70.085 


1.0050 


782.60 


8.3977 


0.30711 


1019.4 


4.5025 


-155.237 


70.060 



of likelihood function. Again the difference at the minimum points is very small and again 
"wrong" minimum is a little preferrable although the total statistics is extreamly high and 
practically unreachable in real experiments. 

Let us return to the previous level of statistics {L = 10^), but change the threshold factor 
/i = 350 instead of 208. Again two minima (Table lV|) difference is statistically unreliable. 

Evidently the more narrow are the resonances, the less is the influence of width depen- 
dence on energy to the form of cross section. Let us set both widthes large — Fi 2 = 100. In 

Table V: Resonance parameters at the minimum points of likelihood function. 



ax 


mix 


Fix 


bx 


m2x 


F2x 


r2x 




0.99946 


782.61 


8.3931 


0.30013 


1019.4 


4.5202 


151.335 


63.799 


1.0051 


782.61 


8.3935 


0.30777 


1019.4 


4.5203 


-153.908 


63.711 
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Figure 7: Fit result of "data points". Cross section model is described by the formula (j94p . 
^jni) = 46.8/(50 - 7). 




\|/2x, degrees 

Figure 8: Likelihood function plot on the second resonance phase V'2x 

order that cross section at the resonance maxima to decrease not so much let us take larger 
amplitudes: a = 10, 6 = 5. Fig. [7] shows the "data points" and fit result. 

Likelihood function plot on the second resonance phase '?/^2x is presented in Fig. [HI This 
time the minimum values are essentially unequal and the "better" minimum has "correct" 
phase. Table lVT] shows the parameters of resonances at these minimum points, confidence 
level of the first minimum is P5o-7(46.844) = 0.318, the second one — P5o_7(64.428) = 
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Table VI: Resonance parameters at the minimum points of likelihood function. 



ax 


mix 


Tlx 


hx 


m2x 


^2x 






10.050 


782.71 


100.61 


5.0111 


1019.4 


100.22 


154.985 


46.844 


10.848 


783.39 


102.03 


5.9767 


1019.1 


99.816 


-173.990 


64.428 





Figure 9: Result of fit of "data points" . Model of the cross section is described by the formula ()94p . 
X^/hd = 52.05/(53 - 7). 



0.0188. 

For completeness let us consider the intermediate case: Fi 2 = 30, a = 3, 6 = 1. Fig. [9] 
demonstrates the "data points" and fit result. 

Again two minima (Table lVnl) are statistically equivalent (the difference of values 
much less than unit). 

After this numerical experiments we can conclude: if resonance width depends on energy, 



Table VII: Resonance parameters at the two likelihood minimum points. 



ax 


mix 


Fix 


bx 


m2x 


^2x 


r2x 


x' 


3.0103 


782.62 


30.128 


0.9986 


1019.5 


30.013 


155.128 


52.054 


3.0745 


782.63 


30.172 


1.1311 


1019.4 


30.019 


-158.262 


51.714 
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Table VIII: Resonances parameters at the two likelihood function minimum points. 



ax 


mix 




hx 


m2x 


^2x 


'^°2x 


Cx 










0.99805 


782.62 


8.3768 


0.30038 


1019.4 


4.5055 


153.925 


3.0128 


1199.7 


100.75 


28.634 


125.969 


1.0535 


782.62 


8.3768 


0.30323 


1019.4 


4.5055 


30.931 


3.1673 


1199.7 


100.75 


-72.318 


125.969 



then the two minimum points of hkehhood function describe not the same cross section 
function of energy. However the difference can be used for cutting off the false minimum 
only under the favourable conditions: high statistics, coverage of wide energy interval with 
both resonances within, and the widths of the resonances must be compatible with the mass 
difference. The final result can be obtained only after comparison the likelihood function 
values at minima: if the difference of levels is much greater than unit, then one can choose 
better set of phases and amplitudes, otherwise one should involve additional considerations 
for the choice of interference phase. 



D. Three resonances 



The case of three resonances with constant widths: 



a{E) 



2m,\a 

E'^ —m'f+ir irrii 



I 2m2be'^2 _, 2m3ce''^3 



with "true" values of parameters 



mi = 782.6, Ti = 8.4, a = 1, 

m2 = 1019.4, T2 = 4.5, b = 0.3, ip2 = 155°, 

mg = 1200, Tg = 100, c = 3, ip3 = 30°, 



(95) 



(96) 



fi = 350, and integrated luminosity equals L = 10®. Fig. [10] demonstrates "experimental" 
data and fit result. 

Likelihood function plot vs the second resonance phase ip2x is shown in Fig.fTTl 
One can see two minimum points on this plot (Table lVlnl) . 

Despite that we could not derive explicit analytical solutions for the case of three res- 
onances it seems that there are at least two equivalent solutions. Let us check whether 
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Figure 11: Likelihood function plot on the phase of the second resonance phase 'ijj2x 

there are some more solutions scanning the space of two parameters: z/'2x and z/'sa;- All local 
minima are presented in Table llXl There are four minimum points with the same values of 
mass, width and likelihood function value. It is quite a surprize that the second resonance 
phase value are different for all points. It means that we should see the four minima at 
the plot of likelihood function, but we have only two of them. In principle it can be. For 
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Table IX: Resonance parameters at the likelihood function local minimum points. 



ax 


mix 




bx 


m2x 


^2x 


^2x 


Cx 


mzx 








1.0535 


782.62 


8.3768 


0.30323 


1019.4 


4.5055 


30.931 


3.1673 


1199.7 


100.75 


-72.319 


125.969 


1.0508 


782.62 


8.3768 


0.23505 


1019.4 


4.5055 


-144.655 


3.0147 


1199.7 


100.75 


-53.143 


125.969 


0.99805 


782.62 


8.3768 


0.30038 


1019.4 


4.5055 


153.925 


3.0128 


1199.7 


100.75 


28.634 


125.969 


0.99545 


782.62 


8.3768 


0.23284 


1019.4 


4.5055 


-21.661 


2.8676 


1199.7 


100.75 


47.811 


125.969 
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-100 100 
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Figure 12: Likelihood function plot on the phase of the second resonance phase V'2a; 

every new minimization run we take as a starting point the final point of the previous mi- 
nimization. Thus the minimization could converge to "bad" local minimum. Let us try to 
get another plot of likelihood function, starting minimization at every point il)2 closer to the 
known "good" minima (Fig. fT^ . 

Now there are all four minimum points on the plot. However the curve is not smooth, 
so probably not at every point of il)2x the global minimum was achieved, although after 
covergence MINUIT executed command IMPROVE, which tries to seach better minimum. 
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Table X: Resonance parameters at the local minima of likelihood function. 



ax 


mix 




bx 


m2x 


^2x 


^2x 


Cx 


mzx 








1.0497 


782.62 


8.3787 


0.30316 


1019.4 


4.5061 


39.985 


3.1635 


1199.7 


100.79 


-64.366 


127.005 


1.0470 


782.62 


8.3790 


0.23583 


1019.4 


4.5058 


-148.442 


3.0063 


1199.7 


100.71 


-45.063 


126.529 


0.99824 


782.62 


8.3792 


0.30038 


1019.4 


4.5057 


153.914 


3.0122 


1199.7 


100.71 


28.687 


126.522 


0.99559 


782.62 


8.3795 


0.23365 


1019.4 


4.5054 


-34.675 


2.8630 


1199.8 


100.65 


47.871 


126.218 



Table XI: Resonance parameters at the minimum points of likelihood function. 



ax 


mix 




bx 


m2x 


^2x 


r2x 




mzx 








11.685 


783.30 


99.681 


10.503 


1019.2 


100.61 


-166.020 


0.75199 


1199.5 


104.28 


157.901 


38.983 


9.9572 


782.86 


100.31 


10.013 


1019.0 


100.10 


152.298 


3.138 


1195.7 


102.44 


25.172 


37.682 



Let us look which set of minima we can obtain if the resonance width depends on energy: 



2m\a 



2m3ce''^3 



(97) 



Again we get the result that for energy dependent resonance width degeneration disap- 
pears (Table |X|). 

Again for narrow resonances this difference is negligible from statistical point of view. 
Let us change the following parameters of two resonances: 

Li = 100, V2 = 100, a = 10, b = 10. (98) 

"Experimental" data and fit result are shown in Fig.[T3l 

On the likelihood function plot vs ip2x (Fig. fT^ one can see the minima, listed in Table lXll 



Table lXin lists the local minimum points found by scanning angles z/'2x and ip^x- 
During this scan five local minima were found. The difference between lowest minimum 
and "highest" one is significant — Ax^ = 4.27. However the difference between global 
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Figure 13: Result of fit the "experimental" points. Process cross section is described by the 
formula x^/no = 37.7/(43 - 11). 




\j/2„ degrees 



Figure 14: Likelihood function plot on the second resonance phase iIj2x 

minimum and closest one is not so big — Ax^ = 0.28. Statistically these two minima are 
almost equivalent. Nevertheless the global minimum has the resonance parameters closer to 
the "true" ones. 
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Table XII: Resonance parameters at the minimum points of likelihood function. 



ax 


mix 




bx 


m2x 


^2x 


i^2x 


Cx 


m-ix 


Tax 


i'lx 




9.9572 


782.86 


100.31 


10.013 


1019.0 


100.10 


152.298 


3.138 


1195.7 


102.45 


25.172 


37.682 


9.769 


782.84 


100.34 


8.8903 


1019.0 


100.08 


138.747 


0.64948 


1195.5 


102.39 


68.199 


37.964 


9.7644 


782.48 


99.576 


8.9027 


1019.2 


99.970 


141.698 


0.62815 


1198.7 


100.53 


77.349 


41.951 


11.881 


783.31 


99.608 


11.810 


1019.2 


100.63 


-153.260 


3.5041 


1199.8 


104.31 


108.19 


39.608 


11.685 


783.30 


99.681 


10.503 


1019.2 


100.61 


-166.026 


0.75197 


1199.5 


104.28 


157.901 


38.983 



VI. CONCLUSION 

As a result of analytical solution of the problem of parameter definition of the two in- 
terfering resonances by experimental data there was demonstrated that for cross section 
parameterization with constant widths there are always two different solutions (for differ- 
ent sets of resonance parameters one gets the same cross section function of energy). If 
the dependence of resonance width on energy is taken into account, then the degeneration 
disappears, but quantitavely two solutions usually differ very little, and this difference is 
determined by many factors. 

For illustration of analytical conclusions a series of numerical experiments was carried 
out. Above conclusions for the case of two resonances are confirmed, although the statistical 
difference of two solutions is not large even if energy dependence of resonance width is taken 
into account. In every particular case this problem should be investigated separately. 

In the case of three resonances for constant widths there occurred already four equivalent 
solutions with the same likelihood function minimum. Analytical solution of this problem 
appeared too hard due to technical difficulties. For one numeric example the system of 
equations was solved, and four different solutions were derived (see appendix [A]) . One can 
guess that the number of different solutions equals 2""^, where n is the number of resonances. 
In case of any number of resonances the degeneration disappears when the dependence 
of resonance width on energy is taken into account, however for narrow resonances the 
statistical difference between different solutions is usually not significant. 

The author is grateful to A. A. Korol for his remarks and recommendations. 
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Appendix A: NUMBER OF SOLUTIONS FOR n RESONANCES 



Let us consider of n interfering resonances: 



fc=i 



mi 



(Al) 



where mfc, are mass and width of k-i\i resonance, < m^+i, and Ak are some complex 
numbers. 

This function is entirely defined by the location and residues of its irregular points, so 
some other function of the form 

2 



E 



A 



kx 



(A2) 



can be equal to the first function over all region of s only if the system of equations is 
satisfied 



A*^ ■ Akx 



y 



A* ■ A, 



y 

^ m] + iTjmj — ml + iV^mk ' 



3 



1, 



,n 



(A3) 



fc=i " J ' J J I' fc=i J 

If we have only two resonances then the system of equations looks like 

A* A I 2^^lml a* a A* A i 2^^lml A* A 

^Ix^lx H G "^lx^2x — ^1^1 H Q ^1^2, 



A* A 2iT2m2 A* A — A* A 2lE2]2}2 A* A 



(A4) 



where G = m\ + iFimi 



iT2m2. 



Let Ai^ = AiZi, Rl = z^zi, A2X = A2Z2Z1: 

Al ■ [{Rlz2 - 1) A2^^ + {Rl - 1) Al] = 0, 
A*2 ■ [{Rl \Z2\^ -l)A2- ■ (Rlz* - 1) Al] = 0. 

Now we can derive Z2 value from the first equation and substitute to the second one: 

[Rl - 1) 



(A5) 



Z2 



I 1 

2mirii?fA^ "rJ' 



{Rl-1). Rl- 



4mi mori To 



0. 



(A6) 



One can see that there are two solutions: the first one is trivial Rl = 1, Z2 = 1, A\x = Ai, 
A2X = A2, and another solution is 



Rl 



1 + 



2imiTiA2 
AiG 



2imiViAl 4mfr2|A2|^ 



GG* 



(AT) 
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If we take the parameters of resonances from the first hne of Table [TTl 
mi = 782.60, Ti = 8.4116, A = 2mia = 1568.3304, 

(A8) 

m2 = 1019.4, T2 = 4.5093, A2 = Im^he^^'' = -555.7337 + 255.8088z, 

then the second solution should correspond to the second line in this Table: Rf = 1.010302, 
Ai^ = AiRi = 1576.3879, Z2 = 0.6557 + 0.77886z, = A2R1Z2 = -563.6463 - 265.09954?, 
a, = |A,|/(2mi) = 1.00715, = |A2.| /(2m2) = 0.3055, ^2x = arg(A2.) = -154.81°. 
Analytical solution matches numerical one within the accuracy defined by rounding errors. 

Now let us carry out similar procedure in case of three resonances with parameters (the 
first row in Table IIXD : 



mi = 782.62, Ti = 8.3768, Ai = 1.0535 ■ 2mi = 1648.98, 

m2 = 1019.4, T2 = 4.5055, A2 = 0.30323 ■ 2m2e^30-93i° = 530.3056 + 317.771i, (A9) 
ms = 1199.7, T-s = 100.75, A3 = 3.1673 ■ 2m3e-^^2.3i9° ^ 2308.1347 - 7240.6307i 

The system of equations: 

(21.4572 - 33.7750i) Rlz2 - (224.715 - 111.758i) Rjz^ = 3271.14 - 145.533z- 

- 3474.402i?2, 



(AlO) 



374.929/2? \z2\ - (100.699 + 7.32447z) Rjz^Zs + (11.5409 + 18.1660z) Rlz*2 = 

= 285.770 + 10.84162, 

48140.6i2? \z3f - (2251.80 - 163.787i) i??z*^2 - (2702.71 - 1344.14z) Rjz^ = 

= 43186.0 + 1507.93i 

where Rf = \zif. We can derive variable Z2 from the first equation: 

' Z2 = (0.653988 + 6.23783i) Z3 - 46.5605 - 73.2892i + 46.9065+67.05i4» ^ 
(14637.5 + 623.355i) Rj \z3f - (1827002 - 90862.5i) Rjz*- 

- (177596 + 97961. 8z) Rfzs + (168317 - 93261. 3z) z*3 + 
+ (163103 + 99669.8i) Z3 = (All) 
= 5321121 - 67.4339Z - 2824793i?? - 
(45646.2 - 13939.2i) Rj \z3f + (114146 + 158750^) - (116606 + 143303z) z* = 

= 43186.0 + 1507.93i 

We got the system of two equations for 23, z^, but both equations are quadratic. Let us 
introduce R^ = {zsf. From the last equation Z3 can be derived as a linear expression of R3, 
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then we substitute and to another equation and derive the only root of -R3: 



^3 



3705.61Rf -13991. 2flf +19810. 0fif-12464.7-Rf +2940.37 
Rf-(i?f-0.892805) ' 
-(290.536+856. 587i)-Rf +(871. 330+2410. 080-Rf-(867.126+2257.000ii!f +286. 440+703.515^ 



(A12) 



i?,J(Rf-0.892805) 

Now we can use expression i?| = 2:3^3 as a final equation for R\: 
{Rf - 0.944882) {Rl - 1) {Rj - 0.89281) {Rj - 0.89747) {Rj - 0.99480) x 



(A13) 







(i?? - 0.94320) + 0.056332 (^2 _ 0.94447) + 0.027832 

This is a polynomial of degree 20 in variable -Ri, or that of degree 10 in variable Rf, so there 
are 20 formal solutions for Ri or 10 different solutions for R^. But there are only five real 
roots Rl> 0: 



Rl 


Z2 


23 


~ 2mi 


J, _ 1^2. 

~ 2m2 


1p2x = 


„ _ 1^3.1 

~ 2m3 


i^sx = 

^^<^*gK(A3=.) 


1 


1 


1 


1.0535 


0.30323 


30.931° 


3.1673 


-72.319° 


0.97205 


(2.6 + 5. 2i) • 10" 


(8.8- 3. 2i) • lOi" 


1.0241 


1.7-10" 


94.5° 


2.9 • 10" 


-92.8° 


0.9449 


0.4937 - 0.6455i 


-0.4810 + 0.8288i 


0.99543 


0.23284 


-21.663° 


2.8677 


47.810° 


0.9473 


-0.5694 + 0.8771i 


-0.1907 + 0.9858i 


0.99803 


0.30038 


153.92° 


3.0128 


28.632° 


0.9974 


-0.7749 - 0.0598i 


0.9014 + 0.3135i 


1.05076 


0.23505 


-144.65° 


3.0147 


-53.14° 



Four of these solutions match with the parameters of resonances in Table IIXI and one is very 
strange (second row). If we substitute the found solutions to the initial system of equations, 
then the four "legal" solutions satisfy the equations within rounding errors, and "illegal" 
second solution does not satisfy neither second equation nor the third one. Obviously this 
false solution corresponds to the case Z2 = Z3 = 00, which should be denied. 

Let us consider the case of three resonances where one of them has infinite width: 



A 



kx 



(A14) 



The system of equations reads: 



0x1 



\A |2 

I^Ol ! 



A: 



2x 



A _|_ \^ ^fcj; 

0^ "T mf -m| +irfe mfe +irimi 



k=l 
2 



k=l 



Al ■ 



m^-ml+ir^.rrn.+irimi 



k=l 



rni-mi+iVkmi:+iT2m2 



(A15) 
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Here we can choose Aq^ = Aq and Aq^ = Aq^, so we get the system of two equations for two 
complex variables Ai^,A2x- Introduce new notations: 

Aix = AiZi, A2x = A2Z1Z2, pi = \zif . (A16) 

It is very hard job to derive the solution for general case. So let us try to solve this problem 
for the following numerical example: 

Ao = 50, 

mi = 728, Ti = 9, = cosl5° -isinl5°, (A17) 
m2 = 1019, 12 = 4, = 39 ■ (cos 155° - i sin 155°) . 

zi,Z2 can be re- written as functions of pi, and for pi we get the algebraic equation of 
12-th degree: 



(Pl 


-1) 


■ (-3.743751 ■ 10^^ 


+ (Pi 


- 1) ■ (9.058522 ■ 10™+ 


+ (P1 


-1) 


■ (1.511509- 10^^ + 


(Pi- 


1) ■ (-1.333590-10^9+ 


+ (P1 


-1) 


■ (-5.218886 ■ 10^3 


+ (Pi 


- 1) ■ (-8.183609-10^^+ 


+ (P1 


-1) 


■ (-2.813253 ■ 10^° 


+ (Pi 


- 1) ■ (1.405807-1036+ 


+ (P1 


-1) 


■ (9.868724 ■ lO^^ + 


(Pi- 


1) - (1.843597-10^3+ 


+ (P1 


-1) 


■ (-8.587479 ■ lO^^ 


+ (Pi 


-i)))))))))))) = o. 



(A18) 



All roots (both real and complex) are located within circle |pi| = 10^"^ One root is trivial: 
Pl = 1, = ^2 = 1. Using Sturm method |2| for the remaining polynomial of 11-th order, 
there was found, that some roots are doubled, and the polynomial 

(pi - if + 2.49 - 10^(pi - 1)^ - 4.29 - 10^^(pi - 1) - 1.07 - 10^^ (A19) 

was a common divisor for the original polynomial and its derivative. Repeating the Sturm 
procedure for the remaining polynomial of 8-th order, it is possible to check that there are 6 
real roots of this equation. Preliminary localization defined exactly one root between every 
two of the following points: 

- 4.9 - lo^ -2.5 - lo^ o, 4.1 - lo^ 8.2 - lo^ 4.295 - 10", 4.297 - io^\ (A20) 



Negative roots should be rejected, because pi can be only positive. Table IXlTTl presents the 
values of positive roots and some additional information. The third row in this table has 
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Table XIII: Real and positive roots of the equation (jAlSp . 







1 i-X \ 


iIja , des. 




\ \ 


iho, des. 


1.000000 


0.9659 + 0.2588i 


1.0000 


-15.00 


-35.346 - 16.482i 


39.000 


-155.00 


1.000413 


0.9702 - 0.2434i 


1.0002 


-14.08 


-35.350 - 4.08 • lO^i 


4.076 • 10^ 


-90.005 


655336.0 


oo 






0)0 






4.292879 • lO" 


1.3863 - 6.552 • lO^i 


6.552 • 10^ 


-90.000 


-35.766 - 15.569i 


39.008 


-156.48 


4.294653 • lO" 


10507 - 6.552 • lO^i 


6.553 • 10^ 


-89.081 


-10541 - 4.075 ■ lO^i 


4.077 • 10^ 


-91.482 



a(s) 



2500 



2499 




Figure 15: Plot of the cross section (IA14P with amplitude parameters from the first row in Ta- 
ble IMIl 

inappropriate solution, because Zi for it goes to infinity. Zi is the ratio of two polynomials 
and for pi = 655336 polynomial in the denominator equals zero. The last two solution are 
quite unexpected because of high value of amplitudes. Original cross section is shown in 
Fig. [151 Within the same interval of s the ratio of cross section of alternative solution and 
the original cross section was evaluated and occurred to be equal to 1 with high accuracy. 
In order to avoid some digital surprises all these calculations were carried out with high 
accuracy of 150 decimal digits, using REDUCE system j^. In order to illustrate the strange 
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0.01 

lm(S,,(s)) 



-0.01 




-0.01 



row 1 



50 
lm(S„(s)) 



row 2 



0.01 



Re(S„(s)) 




-100 



rows 




-50 







Re(S,,(s)) 



Figure 16: Trajectories 6*12 (■§) for different solutions (marked by tlie number of row in Table IXlIIl) 
for the parameter s G (0.36, 1.21) GeV^. Left picture for the solution in the first row, right picture 
is for all the rest solutions. Step for the trajectory plot equals A-^s = 1 MeV. 

two last solutions, in Fig[T6]the trajectories of the complex function 

S,,(s) = ^i^— + (A21) 

s — mi + iTirrii s — m2 + iT2m2 

are presented on the complex plane. 

So we got two solutions for the case of two resonances and four solutions for the case of 
three resonances. It is not enough to choose the rule for the number A^^ of solutions for n 
resonances: it can be Ns = 2"~^, or Ns = 2{n — 1), or something else. 

In order to check whether it is possible to solve the system of equations in every case (at 
least numerically), let us solve a similar problem with four resonances, but choose the most 
simple input data making easier all calculations: 



(A22) 



mi 


= 1, 




= 1, 


Ai = 


1, 


m2 


= 2, 


r2 


= 1, 


A2 = 


1 +i. 


mg 


= 3, 




= 1, 


A^ = 




m4 


= 4, 


r4 


= 1, 


A^ = 


-l + i 
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The system of equations looks like 

'a* (a -U 4 -U A -U l-3» 4 ^ — 247-21i 



150 



15 ' 
207j-25 



1178-1216* 



(A23) 



^0:2 ■ (^^^xl + ^x2 + ^-g^^xS + ^^^2:4) 

/I* 1 3+6i /I _|_ 3+3i A _L /I _L 3-3i /I 
^a;3 ^'1 5~^2:2 "T ^x3 "I ~^x4) — 7Q 

4* ( i±V2i A I 4+8i 4 I 4+4* 4 14 \ 

If we describe the j-th equation in the form 

4 

k=l 

then the solution of every equation can be written as follows 

= ~ " + = Sj ■ i Qj — — ~r~riT' I ' (A25) 



2 5* 



where 



These four solutions together can be considered as a system of linear equations: 

{— k = i 
' j = l,...,n, (A27) 
Gjk, k 7^ j, 

where 

This system can have non-zero solution only if the determinant is equal to zero. If it is, we 
can use the last three equations to express all amplitudes Axj through the amplitude A^i. 

A _ 3 (876+4908i)C3C4+(6615+2205*)C3 + (2744+392i)C4-12250(l+i) ^ a 
^x2 — 5 ■ 4032C2C3C4-5292C2C3-I96OC2C4-54OOC3C4+IIO25 '^2^x1, 

A _ 21 -(4616+7888t)C2C4+21000»C2+3600(i+2)C4-7875(l+2*) ^ a ( A29) 

^x3 — 50 ■ 4032C2C3C4-5292C2C3-I96OC2C4-54OOC3C4+IIO25 '-'3^x1, I J 

A _ 14 (9882-8874t)C2C3+(3*-l)(3500C2-3375C3)-3150(l+30 ^ a 

— 25 ■ 4032C2C3C4-5292C2C3-I96OC2C4-54OOC3C4+IIO25 *-^4^xl- 

If the variables Cj were the predefined constants, then these expressions would be the set of 
infinite number of solutions with arbitrary A^i- But here Cj depend on Axk via the relation 
(]A28|) . And even more, instead of constraint on the determinant of the system (]A27p we 
can use the first equation of the system (]A23[) . which can be presented in the form: 
Rl = F{C„C„C,) = 'J2fl^y, 

(A30) 

^ 4032C2C3C4-5292C2C3-1960C2C4-5400C3C4+11025 

395736C2C3C4-176400C2C3-70560C2C4-490000C2-419040C3C4-165375C3-70560C4+1102500' 
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Table XIV: Results of the search for the solutions in case of four resonances 



Q2 


Qs 


Q4 


Ri 


A^2 




Axi 


min <I> 


-2.000 


-0.1704 


-0.5899 


1.0000 


1.0000 + l.OOOOi 


0.9999 - 1.0002i 


-1.0000 + l.OOOOi 


7.5 • 10-9 


-2.8900 


-4 • 10-8 


0.6880 


1.0492 


0.9378 + 1.2978f 


1.9073 - 0.3002i 


-2.0751 - 2.2721i 


0.009 


-7-10-6 


0.3203 


-0.3965 


1.0459 


1.0731 + 1.3992f 


0.5103 - 2.5446i 


-0.4493 + 1.5012i 


4.2 • 10-"^ 


-2.9180 


0.1103 


0.6893 


1.0536 


0.9359 + 1.3216i 


1.9127 - 0.4163i 


-2.0985 - 2.2274i 


8.1 • 10-8 


0.9510 


-0.0139 


-0.5047 


2.2669 


-1.8270 - 2.5487i 


-1.7284 + 0.0352i 


1.4877 + 0.2031i 


6.5 • lO-'^ 


0.8089 


-5 • 10-8 


0.6846 


2.4065 


-1.8727- 3. 2613i 


-1.4980 - 1.0233i 


-0.6968 + 3.0950i 


0.018 


0.7616 


0.3213 


-0.4026 


1.0472 


1.0754 + 1.4242i 


0.4838 - 2.6127i 


-0.3946 + 1.5174i 


4.6 • 10"^^ 


2.5546 


0.4321 


0.6691 


1.1033 


0.9392 + 1.8158i 


2.5179 - 2.6845i 


-3.2574 - 0.9483i 


1.3 • lO-'^ 



where Ri = \Axi\ 



The problem looks very much complicated. Let us try to solve it using numeric min- 
imization procedure. The free parameters are complex variables C2, C3, C4. Minimized 
function 



[53(F(C2,C3,C4))r 



4 

i=2 



4 



(A31) 



The true solution is found if the minimum value of $ equals zero. Eight different com- 
binations of signs of square roots provide possible eight solutions. Attempt to use the code 
MINUIT {2! failed because of very complicated function profile. Use of 4| brought more 
success. Table IXIVI presents the results of minimization. 

Despite minimization problems this algorithm allows to localize the solutions, exactly 

2"-^ of them, where n is the number of resonances. To improve the amplitudes values and 

make sure that localization is good enough, one can minimize the function 

2 



4 



^*xj GjkAxk - Rj 



k=l 



(A32) 



starting minimization from the found points. There are 7 free paramerters: -Ri = ^{A^i) 
and real and imaginary parts of Ax2, A^s, A^a {^{Arci) = 0). The result of this operation is 
shown in Table IXV[ One can see that for those cases, where all Qj had non-zero values, the 
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Table XV: Improved parameters of solutions. Minimization of ^' started from the approximation 











min ^ 


1.0000 + 0^'^ 


1.0000 + l.OOOOi 


1.0000 - l.OOOOi 


-1.0000 + l.OOOOi 


1.2 • 10-13 


1.0536 + 0i^ 


0.9361 + 1.3216i 


1.9126 - 0.4174i 


-2.0986 - 2.2262i 


5.8 • 10-*^ 


1.0471 + Qi" 


1.0762 + 1.4232i 


0.4831 - 2.6116i 


-0.3947 + 1.517H 


6.2 • 10-6 


1.0536 + Oi" 


0.9360 + 1.3217i 


1.9126 - 0.4172i 


-2.0985 - 2.2262i 


3.9- 10-11 


2.2528 + 0i'' 


-1.8210 - 2.5396i 


-1.7049 + 0.5456i 


1.4949 + 0.2602i 


1.9 • 10-6 


2.3744 + Qi'^ 


-1.5650 - 3.2183i 


-2.0232 - 1.2148i 


-0.6648 + 3.2157i 


1.3 • 10-5 


1.0472 + Oi'^ 


1.0754 + 1.4242i 


0.4838 - 2.6127i 


-0.3946 + 1.5174i 


9.6- 10-15 


1.1033 + 0i'' 


0.9393 + 1.8152i 


2.5170 - 2.6818i 


-3.2560 - 0.9502i 


7.5 • 10-1° 


Additionally found solutions 


2.4858 + Oi*^ 


-1.3958 -4.2950i 


-4.4264 + 0.2797i 


1.1745 + 3. 4443i 


4.8 • 10-1° 


2.3593 + Oi^ 


-1.8330 - 3.4916i 


-2.5821 + 1.8974i 


1.5754 - 0.5896i 


4.9 • 10-^ 



"Improved solution point matched the initial approximation from the Table IXIVI 

''The minimum point moved avay essentially from the approximation in Table IXTVl and matched the solution 
in the fourth row 

■^The minimum point moved avay essentially from the approximation in Table IXTVl and matched the solution 
in the seventh row 

''The minimum point moved avay essentially from the approximation in Table IXTVl Q2 = 0.811, Qs ~ 0.253, 



Qa = 0.676 

'^Reconstructed values Q2 — 0.676, Qz = 0.461, Q4 ~ 0.653 
•^Reconstructed values Q2 = 0.743, Q3 = 0.397, Q4 = -0.210 



improved points practically match the approximate values of A^j. On contrary, for "bad" 
points (rows 2,3 and 6), the improved values of A^j are rather far from approximation in 
Table IXlVl Furthermore the found solutions in rows 2 and 3 match exactly other solutions, so 
the approximations in Table IXIVI were not close to some new solutions. Starting randomly 
from different points, one can find additional two solutions, presented at the bottom of 
Table IXVl 

This exercise shows that the suggested algorithm cannot localize reliably all solutions of 
this problem. But it supports the rule 2"-i for the number of solutions for n resonances 
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problem. 
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